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Abstract 

In this paper we calculate the entropy production of a relativistic binary mixture of inert dilute gases 
using kinetic theory. For this purpose we use the covariant form of Boltzmann's equation which, when 
suitably transformed, yields a formal expression for such quantity. Its physical meaning is extracted when 
the distribution function is expanded in the gradients using the well-known Chapman-Enskog method. 
Retaining the terms to first order, consistently with Linear Irreversible Thermodynamics we show that 
indeed, the entropy production can be expressed as a bilinear form of products between the fluxes and 
their corresponding forces. The implications of this result are thoroughly discussed. 
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I. INTRODUCTION 



In a recent paper [l| we have studied two important thermodynamic aspects of a relativistic binary 
mixture of inert dilute gases using the principles of kinetic theory. The first one concerns the so 
called cross effects, in this case when local thermal equilibrium is assumed, they are the well 
known Dufour and Soret effects The second and most relevant one, concerns with the validity 
of the Onsager Reciprocity Relations (ORR). As we showed in that paper they hold true in two 
representations, or choices, of fluxes and forces. In the first representation, which is referred to 
it in the literature j^l, the heat flux is coupled to a modified Fourier-like force involving both, 
temperature and pressure gradients. In such representation however, the Dufour and Soret effects 
do not appear in their canonical form. The second representation is rather singular. Introducing 
the concept of a "volumetric flow" which arises from the relativistic non-invariance of the volume 
elements in the fluid and whose force turns out to be the pressure gradient, the ORR are shown 
to hold true and the canonical form of such effects is recovered. Further, this representation is 
strictly valid only in the relativistic case. 

As it is well-known in Linear Irreversible Thermodynamics (LIT) Q the appropriate choice of forces 
and fluxes is strongly suggested by the entropy balance equation which in kinetic theory arises 
by simply multiplying Boltzmann's equation by the logarithm of the single particle distribution 
function and averaging over all the velocities of the particles. This procedure allows one to identify 
the entropy production, usually denoted by a as a bilinear form of products between forces and 
fluxes. Symbolically, if Xi is the force associated with a flux Jj, 



where, for an isotropic system, denotes full contraction of the two tensors, necessarily of the 
same rank. 

The full derivation of Eq. (11.11) for a relativistic binary mixture of inert, diluted gases is the main 
objective of this paper. The outcome of this derivation should provide full support for the force-flux 
representations used in Since the whole scheme of both papers is restricted to the tenets of LIT 
we will only need to carry out the calculation using the well-known Chapman-Enskog expansion to 




(LI) 
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first order in the gradients. As shown in Ref. JsJ this procedure is perfectly valid in a relativistic 
framework characterized by Minkowski's metric namely, in special relativity. 

To facilitate the reading of this paper we use the following conventions. Tensors with space and 
time components are labeled with greek subscripts: a,/3 = 1,2,3,4, where the fourth component 
refers to time whereas spatial components are labeled with latin subscripts i,j = 1, 2, 3. Einstein's 
sum convention is adopted for both types of subscripts throughout, but not for sums characterizing 
the species of the mixture. The Minkowski metric has signature {+ + H — }, colons and semicolons 
denote ordinary and covariant derivatives respectively. 

The structure of this paper is as follows: In section 2 the basic ideas of relativistic kinetic theory 
are explained up to a linearized version of the Boltzmann equation consistent with the Chapman- 
Enskog expansion. The solution to such equation is proposed in section 3 using the representation 
introduced in Ref. l|. The general structure of the entropy balance equation and thus the 
expression for the entropy production o are also shown in section 3. Section 4 is devoted to 
establishing the required form of a within this representation (LIT). Lastly, in section 5 we give 
some final remarks. 



II. METHOD: RELATIVISTIC KINETIC THEORY 



The Boltzmann equation |6j for a mixture of two non-reacting dilute species in thermal local 
equilibrium reads, 



'(») 



f(i),a = ^2 J ( 



(ij) 



(HI) 



where the collisional term is given by J], 



j (f(i)f(j) ~ h)fU)) F (ij)^j) dn (ji) dv (jy 



(11.2) 



Here, F^, E(jj) and dTt^^ denote the invariant flux, the invariant differential elastic cross-section 
and the element of solid angle that characterizes a binary collision between particles of the same 
species as well as those between different species. The differential dv*^ stands for rf J[ (>) , which is 
also an invariant. The cross-section E(jj) has special symmetries {?j-9| that guarantee the existence 
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of inverse collisions such that the principle of microscopic reversibility and therefore an H theorem 
are satisfied. 

It is important to notice that the molecular velocity is measured by an observer in an arbitrary 
frame, in which the four- velocity of the fluid is represented by U^. This frame is called the 
laboratory frame. 

From Eq. (III. II) one can obtain the balance equations by multiplying it by the collisional invariants 
namely, the rest mass m^, the four-momentum m^v^ and then integrating over the velocities 
dvZy The complete set of equations can be found in Ref. jlQ|. In this work we only need them 
at lowest order in the gradients. This is accomplished through the use of the Chapman-Enskog 
method of solution to Eq. (III. II) . The particle number density conservation for each species 
reads as, 

n (i) *7£ + U a n^ a = 0, (II.3) 

the momentum balance as 

pU^U^ + h Pv p iV = (II.4) 

and finally, the energy conservation as 

nU u e tU = -pU^. (11.5) 

Here, 

n = ri(i) + ri(j) (II. 6) 

is the total particle density, and 

P = Yl m (i) n (i) G i z (i)) = PW + P{2)> ( IL7 ) 

i 

with 



K 3 

G ( %) ) = — >Tf ( IL8 ) 



IC 2 ( A 

and h^ a = g^ a + c~ 2 U^U a is a projector in the direction orthogonal to U a . Here K, n I — J is the 

\ Z (i) J 

modified Bessel function of the second kind for the integer n. 
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We now expand Eq. fill. II) using the well-known Chapman-Enskog series up to first order 

in the gradients namely, 

/(0 = /(o ( X + *«) > ( IL9 ) 



where is the local equilibrium distribution namely, Jiittner's distribution 



by, 



f 

■J a 



(o) 



n 



(0 



(i) 



exp 



'(00 



12L Il3| which is given 
(11.10) 



with 



k B T 



Moreover, we must now introduce the functional hypothesis namely 



f(i) [% a , v (i) \ n (i), U a , Tj , implying that the representation chosen is defied by the locally conserved 
variables nu\, U" and T. We remind the reader that this assumption constitutes one possibility of 
extracting from the manifold of the possible solutions of the Boltzmann equations those which are 
consistent with the hydrodynamics of the fluid Ref. 14J (also known in the literature as Hilbert's 
paradox). 

From now on, we will develop all the calculations in the local co-moving frame, where the spatial 
components of the hydrodynamical four-velocity vanish, i.e. U m = 0. This frame has the advantage 
that allows us to isolate the purely kinetic effects of the motion of the particles from the convective 



effects 



15Nl7|. Then we can transform all the quantities measured in such frame to an arbitrary 



one with four-velocity U 13 through a Lorentz transformation. Indeed, denoting by the Lorentz 
transformation, the molecular four-velocity in a moving frame reads as, 



v (i) - L -v JX {i)- 



(11.11) 



Here KY^ is the four-velocity in the local co-moving frame. In the classical framework it is precisely 
the definition of the well known peculiar or thermal velocity ll| . 



The definition of the dissipative mass flux, heat flux anc 
obtain the complete set of transport equations (see Ref. 
are given by, 



viscous tensor are established when we 



(0 



10| ) in the local co- moving frame. They 
KftfadKfo, (11.12) 



(i) 



(11.13) 
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and 



7T 



(0 



m (0 | KftKfafadKfo 



(11.14) 



1 - 



-1/2 



respectively. Here 7fc (i) = 

After some algebra, the expansion of the Boltzmann equation ( III. II) with the help of Eqs. ( 111.91) 
( III.3j) . (III.4j) . (1II.5|) leads to the following equation, 



A' 



(0 {-7 %) ^,m + (lnn W ) im + [l + ^ ( 7%) - G (z (0 ))] (lnT),„ 



;n.i5) 



where, 



^(i)^(i)n 



A ( 7 ) A( i)ri 



;n.i6) 



Here T(i) corresponds to the trace of the tensor KJ^.K^ n and it is related to the dynamic pressure, 
while ^A^A(j) n ^ denotes the symmetric traceless part. The collisional linearized kernels are 

[C (0 W ) + C (0 W + <f> U) )] =J2 j ■■■ J f{i)f(j) (0(<) ' + hi) ' ~ ~ ^0')) F (ij)^{ij)dV(ji)dK* U) . 

(11.17) 



The form of Eq. ( III. 151) is crucial in order to identify the thermodynamical forces. As we have 
said, the subtle decision of how it can be rearranged has been studied in a previous work jl]| by 
addressing the necessity of complying with the ORR. This is accomplished through the introduction 
of a pseudo-flux which arises strictly from relativistic considerations and is directly related to the 
term containing the pressure gradient (~ p >m , recall that p^ = n^ksT). Such a flux, we insist 
appears only in the relativistic kinetic theory with p^ m acting as its direct force and can be explained 
by taking the averages of microscopic Lorentz deformations of spatial cells and thus we may refer 
to it as a "volume flux". This idea has never been dealt within the literature before Jj], [jijl and has 
the advantage of allowing a clear definition of both the Soret and Dufour effects in a relativistic 
framework. 
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III. LINEAR THEORY 



By following the arguments that we have discussed in the previous section it is possible to rearrange 
the linearized Boltzmann equation (III. 15j) as follows 



+^ (*w*co») u -m + (ni-i) 

= [<7(0 (i) )+C(<%) + (i) )], 
for species z. Recall that there is a similar equation for species j. Here 

Vm = K(t) m — V (j)m ~ Z [ i - i - LZ ) 

m {j) p P{i) 

represents a new relativistic thermodynamic pseudo-force Vu\ m related with Lorentz contractions 
of the mean free path of the particles. In Ref. Jl]| it has been clearly shown that the corresponding 
transport coefficients satisfy the symmetries required by the Onsager reciprocity relations. Notice 
also that here we have a relativistic generalization of the diffusive force 

A A - rl - r, ( m W G M ~ m U G ( Z U) \ P,m , " / N /ttt on 

d m — a m (i) — ~ cl m(j) — Tiu) 1 {n i0 ) , [LLLS) 

\ p I p n(i) 



with Tiio — — • in the non-relativistic limit one recovers the usual expression 



d m — (mu) - m w ) Vp + Vn i0 . (III.4) 

PP 7 7l(i) 



Now we proceed to the solution of Eq. (IHI.lj) 

7 1 



= -K^A {l) A -f - K%B®V m - K^D®d m - Ljy^. (III.5) 



Notice that when Eq. (1111.51) is substituted into Eqs. (III. 12j) and (III. 13j) we obtain the heat and 



mass fluxes in which the Soret and Dufour effects are clearly identified. Indeed, de Soret coefficient 
is given by, 

J2 f f$KfoK nW A w dKfo (IIL6) 
(0 J 
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and the Dufour by, 



E 

(0,0') 



fti)^; (t* w " G (z W )) K^K^D^dK^ 

z w 



(III.7) 



which have been shown to be symmetric in Ref. 

As it has been thoroughly discussed in Ref. jlj], the pseudo force V m corresponds to a relative flux 
directly associated with Lorentz deformations in the microscopic geometrical aspects of the system. 
Such a flux which we have called a "volume flux" is a novel quantity. Although in non-relativistic 



fluid dynamics the concept of a volume flow dates back to Burnett [20j, it has been recently revived 
by Brenner 21 23|, but in our previous work [l| it arises strictly from relativistic considerations. 



In this paper, such flux which we will denote as Jy is precisely the conjugate of the pseudo force 
V m . Also it is important to underline the fact that both J™ and V m vanish in the non-relativistic 
limit. 



Now we will use Eq. ( IIII.5I) to find the entropy production which, in analogy with the non- 
relativistic case is identified from a balance equation of the form 





Of 



[ns) + V • (J s ) = a. 



(111.8) 



To accomplish this task, we consider the entropy four-flux in an arbitrary frame 
which we decompose as j^l, 



where M is a four vector orthogonal to The invariant a can be thus expressed as 

s»u„ 



or 



(0 

Recalling the properties of invariants, we notice that, 

U^v^ = -7fc w c 2 . 
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(III.9) 



(111.10) 



;m.ii) 



(111.12) 



[111.13) 



Moreover, using that dv%. = d 3 v^)/ , -f v ^ ) = d 3 K^/ / -fk^ i) , we have that 

a = -k B J2 [ hi) (ln/(0 " 1) & K ®, ( IIL14 ) 
(i) 7 

which is readily identified with the local entropy density s for the mixture. In appendix A, it is 
shown that the entropy four-flux defined in Eq. 9 j) satisfies a balance equation of the form 

S^ = a (III. 15) 

where a is the entropy production given by 

a = -k B ^ I J (/(») f(j)) In /(,) . (111.16) 



In the next section we will provide thermodynamical content to Eq. (lIII.16p . Indeed, what we will 



show is that a may be related to Eq. ( 11.11) only when is written in terms of the state variables. 



IV. DERIVATION OF a 



In order to carry out the program outlined in the previous section, we start with the property 
that the Boltzmann equation is a relativistic invariant. We next choose the co-moving frame for 
performing the calculations which will be carried out only to first order in gradients. For this 
purpose we analyze Eq. ( 1III.16|) expanded with the help of Eq. (III.9j) . Thus, 



a 



(0,0") 



,(o) 



Now we expand ln /Jv (l + in a neighborhood of « 1, 



In 



+ ^ln/ ( ^+0 w +O(0 w ) z , 



(o) 



(IV.l) 



(IV.2) 



so we have, 



(T 



(0,(3) 



f (0) f (0) / , 



(7) 



0')) ( m 4) + ^(Oj F (v) E (ij) dn tti) dK (j) dK h- 

(IV.3) 



Noticing that In fj?) is a combination of all collisional invariants, all integrals associated to it 



vanish. Whence 



<y = -k B 

(O.C?) 



^^d^dK^dK^, (IV.4) 



which, with the help of Eq. (111.17p . is easily written as 

a = -k B J2 f [C (^ ( i)) +C(0 (i) 
Finally, by using Eq. (111.150 we get, 



i)dK* {i) . 



(IV.5) 



<7 = -y 



i)dK* i} . 



(IV.6) 



Where is given by Eq. 0111.50 . 

Carrying out the ensuing algebra (see appendix B) one arrives at, 



a 



k 



B 



= ~J m * [d m ] ~ q 7 



T 



T 



— T m V 



K n U -rU. m . 



In Eq. pV.7|) . we identify the vector fluxes namely, the diffusive mass flux as, 



the energy transport, 



where 



W 0) 



Hi) 



Z(i) 



is the specific enthalpy, and the new ingredient, the "volume flux", 



(IV.7) 



(IV.8) 



(IV.9) 



(IV. 10) 



(IV.ll) 



2 

where = P ^' )C depends of the enthalpy through Eq. fill. 71) . The volume flux J^L is defined 
trough a corresponding balance equation |l|. 
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On the other hand we have the tensor fluxes, 



< = = I (IV. 12) 

(i) 

r = E r « = E / ^)^/»V)^(V ( IV - 13 ) 

(0 (*) 

It is important to notice that the structure of these fluxes is the necessary required to obtain the 
canonical form for the entropy production namely, 

a = J2 J i® X u (IV. 14) 

i 

where © is the contraction to a scalar of fluxes with their corresponding forces in accordance with 
Curie's theorem. Equation (1IV.14|) is in complete accordance with the Linear Irreversible Thermo- 
dynamics. Notice that a completely equivalent argument can be given for the other representation 



but since it has been carried out in the literature 



4J we omit it. 



V. FINAL REMARKS 

In this work we have obtained the entropy production for a relativistic binary mixture to first order 
in the gradients using the completely new idea of "volume flux" Ref. jlj|. In such paper, the "volume 
flux" is produced with the next simple idea: Imagine the motion of a single particle, then construct 
an imaginary volume around it with radius a, where a is the mean free path. Such a volume would 
remain spherical in the non-relativistic scheme, but because of the Lorentz contraction it would 
be deformed in the relativistic framework. When we average these microscopic deformations per 
particle we obtain the "volume flux". Indeed, when we multiply Boltzmann's equation by the 
microscopic change in the volume Q7fc (i) , and then integrate over the velocities dKjU we find, 



(/ 



TkMKfofadKfo) = J 7%) (J(u) + %)) dKfo (V.l) 

^vol 1 

which is a balance equation for the change in the volume in the gas, and defines, 

Jvd) = I im K hhi) dK h- (v.2) 
n 



Notice that in the non-relativistic limit, the right hand side vanishes, implying that there is no such 
change in volume. On the other hand, in the one-component limit J" turns out to be a multiple 
of the heat flux q a , indeed t 2 ^ = - J?. 

The corresponding expression for the "volume flux" is valid only in the co-moving frame. This 
however has no restriction since to obtain the same quantity in an arbitrary frame one may simply 
resort to the well-known Lorenz transformations. 

Tfle form of the fluxes We obtained io E q , mM » - D is in 

accordance with LIT and supports the definitions obtained in a previous work Notice that, in 
the non-relativistic limit, the term corresponding to the volume flux Jy in Eq. (1IV.7|) or (IV. 2 II) 
vanishes because 

7fc(0 - G (z (0 ) 0, (V.3) 
so we recover the classical expression for the entropy production 

The generalization of the Soret and Dufour effects found in a previous work Ref. jl| are now 
formally sustained by Eq. ( 1IV.7I) . Their corresponding coefficients are easily found when Eq. 
(IIV.8I) and (IIV.9I) are expanded using the Chapman-Enskog method. Additionally we have two 
more new coefficients related with V m whose physical meaning remains to be studied. 

We would like to emphasize at this stage that the entropy production as defined here and in general, 



in LIT is unfortunate 



25j. Logically speaking it has no meaning since entropy as any other state 



variable such as energy, pressure, volume, etc, can not be "produced". It is a pity that the orig inal 



26| has 



concept of uncompensated heat defined as Ta as originally introduced by Clausius Ref. 
not been kept. Uncompensated heat is the energy that arises in any thermodynamic process due 
to dissipative effects, and further what one can measure in the laboratory is heat, not entropy. 
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Appendix A 



Let 



be the entropy four-flux. Then, 



5 £ = o- 



where 



o- = ^2 J (f(i)fu)) ^f{i) dv h > 



and J {f(i)f(j)) is defined in Eq. (111.2j) . 

Let us start with the four-divergence of Eq. (IV. 4ft . 



Sft. 



kB J2 v (i)h) ( ln ^) ~ l ) dv h 
W 7 



(0 



J ;/i 



« 



dv 



(0 



(0 



(V.4) 



(V.5) 



(V.6) 



(V.7) 



Now we substitute the Boltzmann Equation (III. II) in the second term of the right hand side, which 
reads 



(0 17 



v 



(?) 



^ ( *), 



(V.8) 



then 



so the entropy four-flux is given by 



(*) 



(V.10) 
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and the production term, 



o = -k B J ln f(i) dv h- 

(i),U) 



(v.n) 



It is important to underline here that the identification of S A with the local entropy, S m with the 
entropy diffusive flux and a with the non compensated heat is not complete until the solution 
is determined. Indeed, while remains unknown, none of these quantities can be taken as a 
function of thermodynamic variables because they do not appear in f^y 

On the other hand, 

a = ~ k * E f J (Ao/w) (V-12) 

which, by using the same transformations as those performed in the proof the H theorem it can 
be written as 

ff= \ k *T,f (WW ~ *<> Afl) ln ^^F^a^dn^dv^. (V.14) 



(<).(?) 

Therefore, with the Klein's inequality, 



>0, 



(V.15) 



one obtains 

a > 0. (V.16) 

We recall the reader that Eq. (IV. 16|) is valid for any exact solution of the Boltzmann equation and 
has no inherent physical meaning until one establishes the form of /(j). 

Appendix B 



First recall Eq. flTTO . 



a 



(i) J 



(0) 

(i) 



{dm + ^ ( 7 * (0 - G ^ - ( 7fc(i) - G {Z {i) )) V {i)m ) 

K%yK {i)n )u? m + T ({) UZ 



jm 
;m 



li) dK *i), 



(V.17) 
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so, 

T 



E w I K ™) (ik^ - G ( z d)) ) h) dK h v <i)™ (v. is) 



Here, we identify the bilinear form: fluxes times forces structure. The coefficient of the diffusive 
force is 

I KfttQdKfo. (V.19) 



m(i) 

For the temperature gradient, the coefficient is 



E / K ™)^(lh»- G ^))^ dK h =E W [ / K^-^dK, 

g ho) 



V-fc W -^), (V.20) 



where h(i) = ^y-G ( z (t)) is the specific enthalpy per species. 
The coefficient of the new force V(j) m is given by 

E / K (0 (>) " G ^(O)) ^(0^(0 = E [ / ^(0^(^(0^(0 
(0 (0 

-G (z (i) ) / KftfadKfo 



T^V m {i) c 2 m {i) J 



2 

where = Pw ° , and J„ is defined in Eq. flV.2j) . 
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